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Abstract 

We study some series expansions for the Lambert W function. We show that known 
asymptotic series converge in both real and complex domains. We establish the precise do- 
mains of convergence and other properties of the series, including asymptotic expressions 
for the expansion coefficients. We introduce an new invariant transformation of the series. 
The transformation contains a parameter whose effect on the domain and rate of conver- 
gence is studied theoretically and numerically. We also give alternate representations of 
the expansion coefficients, which imply a number of combinatorial identities. 



1 Introduction 

The equation y a e y = x was solved by de Bruijn and by Comtet [3 [7], as an asymptotic 
expansion: 
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] denotes a Stirling cycle number [HI [TO] . The function u obeys the funda- 



«(^^) = EE(- 1 ) n 

n=l m=l 

and where \ 

L 

mental relation 

1 - e~ u + au - r = . (1.4) 

By rewriting (11.41) in terms of£ = l/(l + cr),a second series for u was obtained by [T7] in terms 
of 2-associated Stirling partition numbers [11] . also called Stirling subset numbers or Stirling 
numbers of the 2nd kind. 
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Since y a e y = x can be solved in terms of the Lambert W function, as y = aW(x 1 ^ a /a), the 
series solutions above are also series for the principal branch of W. We recall that branches are 
denoted W4, and of these only Wo and W-i take real values, obeying the inequalities [8] 

— 1 < Wq(x) < oo for — 1/e < x < oo , 

- oo < W-i(x) < -1 for - 1/e < x < . 

Here, we are concerned only with the principal branch k — 0, although some discussion of series 
expansions for branch k = — 1 is given in section [5j 

The convergence properties of the asymptotic series (11.3p and ( j 1.5ft were first studied in the 
real case in [TTJ, and bounds on the domains of convergence obtained. Here we establish the 
precise domains of convergence of the series both on the real line and in the complex plane. We 
analyze the differences in the properties of the series and find asymptotic expressions for the 
expansion coefficients in (11. 5p . We also consider invariant transformations of the above series. 
The transformations contain a parameter p (related to a above) and retain the basic series 
structure while p varies. The parameter changes the convergence domains of the series as well 
as the rates of convergence. 

A third series is studied, derived from one for the Wright u function [9] . 

m— 1 /; 

m 
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m=l y u/ fc=0 " 

where a = 1/lna; and ujq denotes the Omega constant W(l) = 0.56714329... The notation 
and definition of the second-order Eulerian numbers may be found in [Til [TO] . We give three 
new representations of the expansion coefficients of this series as well as their asymptotic 
estimates. It is also shown that the series (11.51) can also be expressed using the second-order 
Eulerian numbers. Some combinatorial identities follow from the different forms for coefficients, 
including the Carlitz-Riordan identities. 



2 First series expansion (11.31 ) 



It is convenient to consider separately the case in which x > 1, making a, t and u all real. This 
was the case studied in [17], where bounds on the convergence domain were given, and is the 
subject of the next subsection. 

2.1 Convergence for real values 

We begin by stating a convergence criterion in terms of the variables a and r. 
Theorem 1. The domain of convergence of the series (11.3D is defined by the inequality 

lno- < 1 - - + 3£W_i (-e-- 1 ) . (2.1) 



Proof. We rewrite the fundamental relation (II. 4ft by introducing A = rja to play the role of a 
parameter. 

G X ((T, u) = 1 - e~ u + ou - o\ = . (2.2) 
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By the implicit function theorem [21], for fixed A G R, equation (I2.2p determines a function 

MA ( ( x) = ^c m (A)a m (2.3) 

m 

with initial condition u\(0) = in a domain where dG\(a,u)/du = e~ u + <r / 0. The radius 
of convergence of this power series equals the distance from the origin in the complex cx-plane 
to the closest singular point |26|.[3l p. 175, Theorem 4.3.2]. The critical points in the complex 
cx-plane, where dG\(o~,u)/du = 0, satisfy the relations, after using (12. 2p . 

e~ u = -ex , (2.4) 
\-u-l = 1/(7, (2.5) 
-e A -V = e 1/a . (2.6) 

From these, a can be written in terms of the Lambert W function 

<r m = l/WU-e*- 1 ) , (2.7) 

where the mth root is defined by the mth branch of W. Further, by (12.51) . S(l/cx) = $s(— u), and 
for the principal branch 9? (it) G (— 7r,7r). Thus we conclude that there are only two acceptable 
values for m, i.e. m = —1, 0. We now write (I2.6P as 

cr m = -exp{l-A + ^ m (-e A - 1 )} (m = -l,0), (2.8) 

which defines the domain of convergence by the inequality 

lerl < min I — exp {l — A + W m (—e x ~ 1 ) \ I . 

me{-l,0} 1 L " ' 1 

or 

In lerl < 1 - SRA + min UW^-e^ 1 ) . (2.9) 

m6{-l,0} 

Since 3lW-i(x) < KW (x) for all x G R [S], after substituting A = r/cr the condition 02.11) 
follows. □ 

To express the condition ( 12. ip of convergence of the series ( II. 3p in terms of independent 
variable x in ( II. 2p it is convenient to prove the following lemma. 

Lemma 1. Solution of inequality 3?H / _i(x) > a for x < 0, where a is constant, is given by 

I ae a , a < — 1 . . 

x < b = { ' - (2.10) 

I — e a rj csc?7o, a > — 1 

where r] Q G (0, n) is the root of equation r] cot r] Q = —a. 

Proof. We set W_i(x) = £, + if] for real negative x where £ < —1, 77 = for — 1/e < x < and 
£ > —1, — 7r < 1] < for x < —1/e. Then £, r\ obey [8] 



e^(£ cos 77 — 77 sin 77), = ^(77 cos 77 + £ sin 77) . 
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From these equations, one can find the dependence of £ on x explicitly for — 1/e < x < 

i = W_ 1 (x) (2.11) 

and parametrically for x < —1/e 

x = -r/csc^e"^ *' 7 , (2.12) 
£=-770*77, (2.13) 

where — 7r < r/ < 0. 

Now we consider inequality £ > a in two cases comparing a with value —1. 
When a < — 1 the inequality £ > a holds for all x < —1/e because in this case £ > —1 by 
(I2.13p . For —1/e < x < we solve inequality W-\(x) > a due to (12. lip with the result 
— 1/e < x < ae a . Thus £ > a for x < ae a . When a > — 1 the inequality £ > a can have a 
solution only for x < —1/e because £ < — 1 for the rest x. According to (I2.12p and (12.131) the 
solution is given by x < — r] csc(?7o) exp(— r] cot r/o) where r] G (— 7r, 0) satisfies the equation 
—t] cot r] = a due to which the solution can also be written as x < — e a ^ esc r] and r] e (0, 7r). 
Joining both cases, the lemma follows. □ 

Note that in the formula (12.101) . when a > — 1 but a ^ 0, we can also write b = ae a / cos rj . 

Theorem 2. The series (II. 3p is convergent when 

((e/a) a , 0<a<l . . 

x > Xa = { v 1 J ' - 2.14 

gOIJOCSMJO a > J 



where r] satisfies the equation r] cotr] Q — 1 — In a (0 < r] < n), and the series is divergent 
otherwise. 

Proof. We consider the condition of convergence of the series (I1.3P established by Theorem [TJ in 
the real case, i.e. when a > and x > 1. Substituting the expressions f ll.2j) in (12. ip we obtain 

EWLi (-— J > lna- 1 . (2.15) 



Applying Lemma [T] to the inequality (12 . 151) we come to (12.141) . where x a > 1, which justifies 
the assumption x > 1. Thus the theorem is completely proved. □ 

Note. The statement of Theorem [2] was independently reported by A.J.E.M. Janssen and 
J.S.H. van Leeuwaarden [T6] . 

Remark 3. In the formula (I2.14p . when a > 1 but a ^ e, we can also write x a = (e/a) asecvo . 

Remark 4. It follows from (12.141) that x a — > 1 and x a — > oo as a — > oo. In addition, 

one can show that a function of a defined by (12. 14|) zs monotone increasing. Therefore, the 
larger a, the less the domain of convergence of the series (1 1 . 3 j) . 
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Corollary 5. Series (jl.3p for a = 1, i.e. for W function, is convergent when 



5R 



W-i 



\nx 



>-l , 



which is equivalent to x > e, and divergent when the opposite inequalities hold. 
Proof. Follows immediately from (I2.15P and (12.141) . 

Remark 6. In terms of variable a = 1/lnx series ( 11. 3ft is convergent for < a < 1. 



(2.16) 



□ 



We now prove a statement, relating to divergence of the series ( 11 .3j) . which was found by us 
earlier than the conditions ( 12. 14ft but unlike Theorem |2] concerning positive a it deals with any 
a^O. In addition, the statement demonstrates an interesting application of the ratio test to 
the series (II. 3p . 



Theorem 7. The series ( 11. 3D is divergent at least for 



e~ H <x< e bH 



(2.17) 



where 



W (1/ \a\) when \a\ < 1/e 
1 when \a\ > 1/e 



Proof. Changing indices for summing the expansion (11. 3p can be written through a double 
series [T7] 

oo oo 

« = ^^c m , z , (2.18) 



where 



c m ,z = c mt i(a, t) 



m=l 1=0 

(-1Y 
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For the column-series c m ^ the ratio test gives 
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m— ¥oo 



c m+l,l 
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m— >od 





'l + m + 1 






i + 1 




(m + 1) 


/ + m 
J + 



as according to pQ 



lim 

p— f oo 



p + 1 
l + l 



p 



V 

l + l 



for fixed I 



5 



in our notations. 

For the row-series ^ c m ^ we have 



lim 

l— too 



c m,l+l 



| er | lim 

l— >oo 



1 + m + l 




1 + 2 






I + m 











because by [TJ 



lim 



/ + m 
/ + 1 



for fixed m. 



l^oo (I + l)2m-2 2 m " 1 (m-l)! 

According to [201 Theorem 2.7] the series (I2.18P (and therefore (11.31) ) is divergent when |r| > 1 
or \a\ > 1. Expressing these inequalities in terms of x by fll.2p and uniting the obtained sets 
we come to the stated inequality (I2.17p . □ 

By Theorem [7] for a = 1 the series (11.31) is divergent at least for e _1 < x < e, which is 
consistent with Corollary [5j 

For comparison, the curves described by equations (12. 14j) and (I2.17P are shown in Figure [1] 
by solid and dash-dotted lines respectively, together with the bound given in [17] . 




Figure 1: Boundary of domain of convergence of series (II .3p given by (I2.14p is depicted by solid 
line; dashed and dash-dotted lines are described by equations given in reference [T7] and (I2.17P 
respectively. 



2.2 Convergence in complex case 

In this subsection we consider equations (ll.ip . (I1.2p only for a = 1 (under the same relation 
(ll.4p ) and derive convergence conditions for the series (11.31) in the complex case using the results 
obtained in Section 12.11 in the real case and providing a continuous continuation of the latter. 

We set 

a = l/\az,T = lnlnz/mz , (2.19) 
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where z = x + iy is a complex variable and In z denotes the principal branch of the natural 
logarithm. Then the right-hand side of the series ( 11. 3p represents a function of the complex 
variable z and the following theorem holds. 

Theorem 8. The domain of convergence of the series ( II. 3p in the complex z -plane is defined by 

mv m > -!> ( 2 -20) 

where the branch W m is chosen as follows 



m 



•1, — 7r < arg z < 
1, < arg z < 7r 



Proof. Repeating the proof of Theorem [T] under assumption A G C we come to an equation 
which is different from (12.91) only by that m G Z. Then substituting (12 . 19[) in there we obtain 
(cf. Q2HBP) 

miri mV m ( > -1 . (2.21) 



tilts 



Now we cut the complex z-plane along the negative real axis and set arg z G (— 7r, it]. We 
consider inequality ( 12 . 2 1 [) in domain D = {z G C| — n < arg z < 0} and assume that there 
exists some value m = q such that the domain of convergence in D is defined by equation 

mV q (-\nz/e) > -1 (2.22) 

and its continuous boundary Jzf is given by 

mV q (-\nz/e) = -1 . (2.23) 

The domain of convergence found in real case is defined in a similar way. Specifically, in domain 
{z G M| 2 > 0} we have q = —1 and the boundary z = e (see Corollary [5]) . We require that 
in the limiting case arg z — > 0— equation ( 12.22ft become equation ( 12.16P and show that there 
is an unique value q = — 1 satisfying this requirement. (If there were several such values of q, 
it would mean that the boundary Jf is composed of several pieces of different curves, and to 
identify them one should reduce domain D, i.e. consider its subdomains.) Separating the real 
and imaginary parts of W q (—hxz/e) = — 1 + ir\ in the usual way, we find 

sin rj — 7] cos 7] = arg z (2.24a) 

cos 77 — (— 77 sin 11 = In I z I (2.24b) 

These equations describe a set of the boundary points which can be found in the following way. 
Given a value for argz one can find rj from (I2.24al) which being substituted in (12.24bp yields 
the corresponding value of ln|z|. However, for fixed arg z G (— ir, 0] the equation ( 12.24a[) has 
an infinite number of solutions. We select a solution to provide a continuous transition to the 
real case when arg z — > 0— and when the boundary of the domain of convergence is defined by 
3?W / _i(— lnz/e) = —1 (cf. (I2.16P ). An elementary analysis of the equation (I2.24al) shows that 
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to meet these requirements one needs to choose a solution of this equation from the interval 
r] G (— 7r, 0] and set q = — 1 in (12.231) . Since by (I2.24al) such solution exists if and only if z G D, 
the above assumption is approved and the domain of convergence in D is described by (I2.22p 
with q = — 1, i.e. 

mV-xi-lnz/e) > -1 . 

Due to the near conjugate symmetry property of W function [8], i.e. Wk(z) = W-k{z) when 
z is not on the branch cut, we obtain the convergence condition 3tWi(— In z/e) > — 1 in the 
domain {z e C| < argz < 7r}. Thus the theorem is completely proved. □ 

Remark 9. The 'branch splitting' in the proved formula (I2.20p is due to the branch choices 
for the Lambert W function and similar to the effect that occurs in the series for W about the 
branch point ITU, Sec. 3}. 

Remark 10. The inequality opposite to (12.201) defines the domain where the series (11.31) is 
divergent. This domain is finite (it encloses the origin z = 0) and contains a subdomain 
defined by inequality \o~\ > 1. Therefore, unlike the real case (see Corollary^) in the complex 
case the condition \o~\ < 1 is only necessary but not sufficient for convergence of the series (11.31) . 



3 Series (11.51) 

3.1 Convergence in real case 

We regard the expansion fll.5p as a power series around r = where variable a plays a role of 
a parameter. 

Theorem 11. For a > and a > 0, the radius of convergence of the series (11.51) is exactly 

r*(cr) = |1 + a — alna ± mo\ . (3-1) 
which is equiavalent to the condition of convergence of the series (11.51) as 

|(j(ln a - In a) | < ^(1 + a - a lna) 2 + vrV 2 . (3.2) 
Proof. We rewrite the fundamental relation (11.41) in the form of equation 

F a (r,u) = Q, (3.3) 

where 

F a (t, u) = 1 - e~ u + au - r , (3.4) 

and analyse this equation similarly to that in the proof of Theorem [TJ By Implicit Function 
Theorem the equation (13. 3p determines a function u a (r) = J2m c m( CT ) Tm with initial condi- 
tion w^O) = in a domain where dF a (T, u)/du = e~ u + cr^0. The initial condition is justified 
by dF a (0, 0)/du = 1 + cr ^ 0. Since the critical points are defined by the same equation as in 
Theorem (TJ they are given by (12 .4p and the corresponding values of r are 

rf } = 1 - e- u * ] + auf ] = l + a~a\na + ina(2k - 1), k E Z (3.5) 
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The radius of convergence is equal to the distance from the origin in the complex r-plane 
to the closest singular point [3j Theorem 4.3.2]. Among the critical points ( 13. 5 j) there are two 
the nearest to the origin equidistant points which correspond to k = and k = 1: 

rf°) = 1 + cr — alna — ina, (3.6a) 

T} 1 ' = 1 + cr — crlncr + iixo. (3.6b) 
The corresponding values of are 

wi 0) = - In a - in, (3.7a) 

«f ] = -lna + ivr. (3.7b) 

Since the expansion coefficients of the series ( 11.51) are real, the closest singularities can 
appear as a conjugate pair only [15]. Based on the Weierstrass's preparation theorem [2TJ E] 
we will show that the points ( 13. 6p are singular, each corresponding to a square-root branch 
point of function u = u a (r) in the complex r-plane. We will also find a behavior of function 
u = u a {r) near the points (13. 6p used then for a study of an asymptotic behaviour of the 
expansion coefficients of the series ( 11. 51) . 

Let us consider, for example, point r = T*f\ Expanding the left-hand side of equation (13.31) 
into a Taylor series near the point S (rf \ j we obtain 

F . (S) + d -^l(r- rP) + „<<») + a ^< S) ^ " T ' <0)N 2 



dr v * y v * y <9r 2 2 

drdu [ ' ,[ • I + &u? 2 + 

where dots denote the skipped terms of the higher order. Since 

n dF a (S) d 2 F a (S) , m , 

F *(S) = 0, = 0, —gjr- = - exp {-u[ , and — = -1 

the last equation becomes 

( (o) x 2 
[u — Ui ' 
(n\\ i mn \ 

r - t 



(°>) - exp V ' +.- = . 



Thus, in accordance with the Weierstrass's preparation theorem [211 P-UlL equation (13. 3p is 
locally equivalent to the equation 

(ri°) - r) ~ exp (-„}») 
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It follows that at r = r* function u = u a {r) has a singularity corresponding to a square- root 
branch point as near this point 

u ~ ± exp 

or substituting (I3.7aj) 




2ri 0) / r\* 



a 



u ~ — In a — Z7r ± i\ (l j-- . (3. 



r (o) 



It is not difficult to show that if we consider the values of the function (13. Sp in the interior of 
the circle of radius (13. 1 p remaining in the vicinity of r = ri°^ then the function (13. Sp taken with 
the plus sign only satisfies the condition — tt < < 0, which corresponds to = — tt < 

at point r = ri°^ itself by ( 13 . 7aj) . Moreover, since in the mentioned vicinity —tt < Qt/ct < 0, 
we have — tt < ^sW < tt, which corresponds to the principal branch of W function [8]. Thus we 
come to conclusion that the function u = u a (r) behaves near the singularity (I3.6a|) like 



2ri 0) /. r 



a 



u ~ —In a — iix + i\j — ^— ( 1 ^ ) as T->rf. (3.9a) 



One can show in a similar way that near the singularity ( I3.6bl) the function u = u a (T) behaves 
like 

I 2t W ( T V 

u \na + i7r-i\ — — 1 m ) as r . (3.9b) 



a 



' * 



Thus the points ( 13.61) are singular and we immediately obtain expression (13.11) ; the inequality 
( 13. 2 p follows from ( 13. ip as r = — cr(lncr — In a) due to (1 1 . 2 [) . The theorem is completely 
proved. □ 

Remark 12. From the values ( 13. 6 p and ( 13. 7ft . we /inc? Vl/(x) = — 1 /or &ot/i A; = and k = 1. 
Although it is well-known that this value of the Lambert W function corresponds to its branch 
point and asymptotics ( 13. 9 p can be obtained immediately from the results in |2l HE/, we derived 
these asymptotic formulae to demonstrate a method based on the Weierstrass's preparation 
theorem. 

The inequality (13. 2p can be written in the form |lna — \no~\ < g(o~), where 



o(a) = v / vr 2 + (l + l/o--lna) 2 , (3.10) 
and solved with respect to a. Then the following theorem follows. 



Theorem 13. Let a c = 1.059945... is the unique root of equation g(a)^(l + 1/cr) 2 — 1/(1 + 
1/cr) = 7r and a c = a c exp(g(a c )) = 41.349171..., where function g (a) is defined by ( 13.101) . Then 
the domain of convergence of series ( 11.51) depending on a is defined as follows. 
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(i) For < a < e, the series (11.51) is convergent when < a < o~ a , where a a is the only root 
of the equation 

In cr — g(a) = In a, (3-11) 

which is equivalent to 

x > x a = e a/aa . (3.12) 
The series is divergent when o > o~ a or 1 < x < x a . 

(ii) For e < a < a c , the series (II. 5p is convergent for any a > or for any x > 1. 

(iii) For a > a c , the series (ll.5p is convergent when o < \i a or a > v a , where /i Q and v a are 
roots of equation \lna — lna| = g(o~), which is equivalent to x > e a ^ a or x < e a ^ Ua . The 
series is divergent when fx a < o < u a or e a ' Va < x < e a ^ a . 

Proof. We give details only for the first part (i). In a particular case a < a, equation (13. ip 
can be written as (13. lip with g(o~) defined by (13. 10p . The left-hand side of equation (13. lip , 
being a monotone increasing function for positive a, goes to — oo and 1 when a tends to 
and oo respectively. Therefore, for < a < e the equation has the unique solution. Applying 
Theorem [TT] the theorem follows. □ 

Corollary 14. The series (II. 5p for W function is convergent for < a < o\, where o\ = 
224.790951... is the only root of the equation 



1-oTnoi = y/(l + a-o-\na) 2 + ii 2 a 2 , (3.13) 

which is equivalent to 

x> Xl = e 1/CTl = 1.004458... . (3.14) 
The series (jl.5p is divergent for o > 0\ or for 1 < x < x±. 



Proof. Follows from Theorem [TBT i) for a — I. □ 

Remark 15. In terms of the variable x the series (II. 5p for W function is convergent for 
x > Xi > 1 rather than for x > 1 though x± is very close to unit. 

Remark 16. Elementary analysis of equation (13. lip with substituting a = a/lnx(a) therein 
shows that x(a) > 1 for < a < e and has one maximum e e " = 1.044161... at point a = 
e _7r /W(l/e) = 0.155186.... This means the dependence x(a) to be a very weak and x a = 
x(a) can be evaluated with a good precision (with the relative error less than 5%) by a simple 
approximate equality x a ~ x\ (0 < a < e) which becomes accurate when a = 1. 

Remark 17. The solution o = 0\ of the equation (I3.13P is much more than unit and can 
be found approximately with a good precision. Specifically, taking square of the both sides of 
(I3.13P and leaving the main terms we obtain a 2 — 2a 2 In a + 7i 2 a 2 — 2a In a ~ 0. Searching for a 
solution of the approximate equation in the form a = exp(i4^-)(l + 5), where the exponential 
factor is an exact solution of the approximate equation with neglected last term and a correction 
term 5 is to be determined, we obtain an approximate value in deficit o~\ ~ exp(i^-) — = 
223.8126969.... Taking into consideration of the terms of higher powers in S in a similar way, 
one can obtain a more accurate value. 
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Remark 18. For fixed a > and a — 1, the singular points (13. 6p of function n^r) correspond 
to those of the Wright function u = u>(s), which are = £ ± m, where £ < —1 (see subsection 
\6.2\) . Indeed, we have r* = 1 — o In a — crs*, z.e. r* = 1 — a In a — £cr =p mcr. Since 3ftr* /jas 
i/je minimum at £ = —1, t/ie closest singular points are defined by equations, which are exactly 

d2D. 

Remark 19. TTie convergence condition (13. 2p /ias a c/ear geometrical interpretation in (cr, r)- 
plane. For example, for a = 1 one can s/iou> £/ia£ m accordance with the inequality ( 13. 2[) . u>/ien 
a < o"i i/ie curve L described by r = —o In a zs located inside the region S bounded by curves 
t = ±a/ (1 + a — a In cr) 2 + 7r 2 cr 2 , which expresses the condition of convergence of the series 
(11. 5p . However, at point o = o\ the curve L leaves the region S through the lower boundary 
curve that can be described for large a by the asymptotic expression 

r(a) = -y/CL + a -a\na) 2 + vrV 2 = -a In a + a - 1 + — - - +1 + ( 

2 In a \ In a 

It follows that afterwards the curve L remains below the lower boundary of S, which corresponds 
to the divergence of the series (II ,5p for o~ > o~±. 

Now we consider case a < 0, which should be done carefully as by Implicit Function Theorem 
it should be dF a (0,0)/du ^ due to the initial condition u a (0) = and therefore the value 
a — — 1 should be excluded. It follows from (I3.5P that when o < and cr ^ -1, i.e. a = \a\ e l7T 
and |o"| 7^ 1 there is only one the nearest to the origin singularity given by f!3.6bp 

7-W = l + o- - a lnH (3.15) 

that lies on the positive real axis. Correspondingly the radius of convergence instead of (13. ip 
is the modulus of the right-hand side of (I3.15p . 

Finally, when a = the series following from (ll.4p 

oo m 

u = -ln(l-r) = V — (3.16) 
z — ' m 

m=l 

is convergent for |r| < 1. 

Note. When a = —1, r* = by (13.151) . i.e. the series diverges everywhere. We also note 
that in all cases considered above the condition of convergence of the series (I1.5P is described 
in an uniform manner, particularly, the radius of convergence is given by (13.11) . 



3.2 Comparison with series (11.31 ) 

Let us compare the domain of convergence for the series (11.31) and (II .5p . Both can be represented 
in the form 

oo 

u=Y,c m {o-)r m (3.17) 

771=1 

(see (I2.18P for the series (II. 3p ). However, by Corollary [TH and Corollary [5] the series (II. 5p has 
a much wider domain of convergence than the series (11.31) (not only in the real case but also 
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in the complex case, see Figure |2] below). To undestand this phenomenon we note that both 
series (11 .3p and (II. 5p are defined for x > 1, therefore the closer the boundary of a domain 
of convergence to 1, the wider domain of convergence. However, the expansion coefficients 
c m (a) in the series ( 11.3)) are given by power series near a = whereas in the series ( 11.5)) . the 
expansion coefficients are defined through function ( = ((a), i.e. c m (o) = c m (C(o)), where 
£(o) = l/(l + cr), and are given by power series near ( = 0. Since a = 1/ lux becomes larger as 
x is approaching 1, the series (jl.5p has a wider domain of convergence. This corresponds to the 
fact that the function ( = ((a) maps the interior of the unit circle \a\ = 1 into an unbounded 
domain which is the right half-plane K.£ > 1/2. We also note that the series ( 11.3)1 and ( 11.5)1 
have common values in the domain where they are both convergent, therefore the series ( 11. 5p 
is an analytic continuation of the series ( II. 3p . 

In terms of variable ( the series (11.5)) becomes [17] 



and can be regarded as a result of applying the Euler's transformation for improvement of 
convergence of series [12]. Indeed, the standard Euler's transformation associated with changing 
variable to extend a domain of convergence of the series ( 11.31) is p = o/(l + o) [22]. Since in 
terms of a new varibale the fundamental relation ( 11 .4p is written as 



it would be natural to introduce variable £ = 1 — p = 1/(1 + 0") rather than p. The series 
( E3$ , @ng) were first found in [T7]. 

One can also show that a representation of W function through the function u T (a) = 
Ylm=i c m{T~)<J m , where r plays a role of parameter, can not extend the domain of convergence 
established for series (1 1 . 5 p . Indeed, in this case equation (I3.3P changes to F T (a, u) = where 
F T (a, u) is still defined by the right-hand side of (13 .4p but with initial condition u T (0) = — ln(l — 
r). By the Implicit Function Theorem it should be dF T (0, — ln(l — r))/du ^ 0, which gives 
r ^ 1, i.e. \t\ < 1, and substituting r = — olno yields < o < 1/uo as a necessary condition 
for convergence (cf. < o < o\ in Corollary (HJ) . 

Thus among the series with the considered structures the series ( 11. 5ft has as wide as possible 
domain of convergence. 

3.3 Asymptotics of expansion coefficients 

Once the behavior of function u = u a (r) near the nearest to the origin singularities has been 
established one can find an asymptotic formula for the expansion coefficients of the series (11.51) 
using the Darboux's theorem about expansions at algebraic singularities [7J 0]. The similar 
approach, based on the Weierstrass's preparation theorem and the Darboux's theorem, was 
applied to asymptotic enumeration of trees in [23] . 

According to the Darboux's theorem and found estimates (13.91) for o > the expansion 




(3.18) 



1-p 



e u + u + r — 1 
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coefficients in the series (13.171) have an asymptotic formula for large m as 




2tI 



(i) 



" r(4)(r«»» 



1 

777,2 



° r(-J)(r. 



(l)\ m 1 

> 777,2 



+ 



3 

WJ2 



(i) 



or 



Cm. W) 



V2 



as 



r (-§) = -2v^F. Setting r, 



(i) 



Tiaras 



r, 



(0)N m ~2 



7* 



(3.19) 



Cm. C 




e* 6 * 1 we find 
2 sin (m — ~) #i 



7TO" 



m-i 3 ' 

r* ma 



as m — > oo 



(3.20) 



where = t*(<t) is defined by (13.11) and Q\ = arg(l + a — a In a + mcr). 
Specifically, for a > 



arctan 



7T 



7T + arctan 



l-lno- + l/(7 5 

7T 



1 - In 



<7 



1/a' 



if < a < 



if a > 



1 



W(l/e)' 



W(l/e)' 



It follows from (I3.20|) that for large m the expansion coefficients in the series (jl.5p disclose 
their oscillatory behavior due to sin function though the amplitude decays as r*(c) > 1 for any 
cr > 0. Since the series (11. 5ft can be interpreted as a result of applying the Euler's transformation 
to the series ( 11. 3D (cf. ( 13.181) ). we note that some cases of oscillatory coefficients resulting from 
the Euler's transformation are studied in [14|. 

In order to find an asymptotic formula in case when a < 0, suffice it to take in (I3.19P only 
the first term with (13.151) 



Cm. C 



a/27t \a\mi (1 — \a\ + \a\ In \a\) m 2 
Finally, for case o = 0, it follows from ( 13 . 1 6[) that for any m G N 



as m oo 



(3.21) 



Cm(0) 



(3.22) 



3.4 Convergence in complex case 
Theorem [TT] is extended to the complex case. 

Theorem 20. For complex a, the radius of convergence of the series (II. 5p for a = 1 is 

|1 + a — erlner — ma\ when < , (3.23a) 



tact) 
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7"*((t) = |1 + a — a In a + nrer| u>/ien 3a > . (3.23b) 

In the complex z-plane this is equiavalent to the series (11.51) for W function is convergent 
everywhere in the exterior of the boundary line defined by equation 

| — cr lncr| = |1 + a — oTna ± ma\ , (3.24) 

where a = 1/ In z and sign minus or plus is taken respectively in the upper or lower half-plane. 

Proof. Repeating the proof of Theorem [11] under assumption a G C we obtain the same equa- 
tions (12.41) and (13. 5p for singular points u* and ri^ respectively, where fceZ. However, many 
of the singular points do not correspond to the principal branch of W function and relate to 
the other branches. We are going to find acceptable values of k for which singular points relate 
to the principal branch of W. 

To find acceptable values of k, we utilize a relation r = — cr In a following from (ll.2p to 
obtain 

u = W(e s ) - 1/cr-lncr , (3.25) 

where s = — In a. Let us consider values of u in the e- vicinity of the point u* . Comparing 
between (I3.25j) and (I2.4p gives 

m{2k - 1) + ee* = W(e s ) - l/a , 

where — tt < tp < it. Setting z = \z\e ld (— n < 9 < n) in cr = 1/ In z and separating the 
imaginary part in the last equation we obtain 

n(2k - 1) + e sin <p = %W - 9 . (3.26) 

Since for the principal branch — 7r < < tt, we find —1/2 — e sin (p/(27r) < k < 3/2 — 

e sin ip/(2n), i.e. acceptable values are k = and k — 1. 

Now we note that both points ri°^ and ri 1 ^ are singular, particularly, they correspond to a 
square- root branch point of function u = u a (r) for the same reason as in the real case (see proof 
of Theorem [TT!) . Taking into account this result we consider equation (I3.26P for e = in two 
cases k = and k = 1. When k = 0, we have = 9 — n. Since — tt < < 7T, only positive 
9 satisfy this equation, i.e. < 9 < n. Similarly when k = 1 we have = 9 + tt which 

holds for — 7r < 9 < 0. Thus we conclude that the curve I— oTn cx| 



is located in the 



upper z-half-plane and the curve |— a In a \ 



is located in the lower ^-half-plane, these 



curves being symmetric with respect to the real axis. Hence, the equation (I3.24p describes the 
boundary of domain of convergence of the series (11.51) in the complex case. In addition, since 
a = (In 1 2 1 — iff) 1 1 In z\ 2 , 9 and are of opposite signs and the equations (13.231) follow. The 
theorem is completely proved. □ 

The curve defined by equation (13.241) in the complex z-plane is depicted in Figure |2] by solid 
line in the upper half-plane only (corresponding to the negative sign) as it is symmetric with 
respect to the real axis. The exterior of this boundary line can be regarded as the domain of 
analytic continuation of the series (II. 5p from the part of the real axis x > x\ (see Corollary 
[T4l to the complex z-plane. For comparison in the same figure it is shown (by dashed line) the 
boundary line of the domain of convergence of the series (11.31) defined by equation (12.201) . 
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Figure 2: The domains of convergence of series (11.51) and (11.31) in the complex z-plane are 
located in the exterior of the curves depicted by solid and dashed lines respectively. 

Remark 21. We note that the case \a\ < 1 reveals a connection between the series (I1.5P and 
(II. 3p . In particular, the case permits to expand 1/(1 + a) in powers of a in the former that 
after some rearrangments can be reduced to the latter [1 7/. In accordance with Theorem 
the series (11.31) is convergent in domain V in the complex a -plane defined by (12.201) (written 
in terms of a = l/\nz). One can show that the domain V is contained in the unit disc 
U = {cr g C | \cx\ < 1} (cf. Remark [70j) with the boundaries of V and U having one common 
point a = 1 (where both series are convergent). The series (11.51) is also convergent in V but 
has a wider domain of convergence being convergent (for \a\ < 1) in U D H where a domain H 
bounded by curve (13.241) . 

In the end of this subsection we give asymptotics for the expansion coefficients of the series 
(ll.5p as m — > oo when ^ 0. It follows from the proof of Theorem [20] that in this case there 
is only one singularity when Ser < and r = 7* when Scr > 0. Therefore, one can 

use formula (13.191) keeping only one corresponding term (unlike case of real a when there occur 
two singularities and both terms constitute the asymptotic formula (I3.20p ). Thus, taking (13.61) 
we find 

. 1 ±i 

CmiC) ~ . ; 7777 aS 171 — V CX3, 

v/2WV 2 (1 + a - a\na ± maf 1 ^ 2 
where sign "+" ("— ") is taken in case of positive (negative) Qa. 



16 



3.5 Representation in terms of Eulerian numbers 

The expansion coefficients of the series (11. 5p can be expressed in terms of the second-order 
Eulerian numbers [Hi [10] . To show that we combine 



u = W(e s 



1 - T 



a 



where 



s = s(a, t) = lna . 

a 



and ( 13.171) . then the coefficients c m (cr) in the right-hand side of ( 13.171) are 



1 d m 
c m {a) = -— -W(e° 
ml dr m 



T=0 



l\ m (f 



ml \ o I ds 



-W(e s 



(3.27) 
(3.28) 

(3.29) 



s=— In (7+1/ a 



as ds/dr = -1/a by ( 13~28|) . 

Because of (I3.27P the formula (I3.29j) is valid for m > 2, for m = 1 we have 



c x {a) = - + ±W (e s 
a dr 



(3.30) 



T=0 



Since ED] 



d r ' 



-W(e s 



q m {W{e s )) 

,2m-l ' 



ds m v ~ ' (\ + w{e s )Y 
where the polynomials q n {f) can be expressed in terms of the second-order Eulerian numbers 

mi ma 



q m [r) 



m— 1 

E 

k=0 



m — 1 



(-l)*r 



and 



we finally obtain 



Ci((7) 



VT(e s 



U= — lncr+l/tr 



W ( -e 1 ^ 
a 



1 

a 



m^.m—1 



1 + 0- 



, c m (cr) 



-1)"V 



m!(l + cr 



i2m- 



m—l 



fc=0 



m — l\\ (— 1 



k 



a 



k+l 



m > 2. 



(3.31) 



Substituting (13.311) into the right-hand side of (13. 17f) results in a desireable formula 
1 + a + 5^ m!(l + (r) 2m 



m — 1 



A- 



m—k ^m—k—2 



(3.32) 



m=2 v ' k=0 

By introducing the variable ( = 1/(1 + cr) the series f!3.32j) can also be written as 



OO yy, m—1 II 

T \ m — 1 



U — T' 



m=2 k=0 



^yn+k^m+k+l^ Q\ 



m—k— 2 



(3.33) 
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We note that the expansion (I3.33|) does not contain terms of the second order in (. 

The series (Il.5p . fl3.18p .f l3.32p . and (I3.33j) have the same properties including the domain of 
convergence and the asymtotic estimates for the expansion coefficients studied in Section 13.31 
This fact leads to some combinatorial consequences considered in Section [7J 



4 Transformed series 

4.1 An invariant transformation 

The above studied expansions (11. 3p and (11.51) are limited in their domain of applicability by 
the fact that a and r are each singular at z = 1, restricting their utility to z > 1. In addition 
to the domain of validity of the variables, there is the question of the domain of convergence 
of the series ascertained in theorems El [HJ [11] and [20] 

In this section we consider transformations of the series (I1.3P and (11. 5p referred as to trans- 
formed series. Our aims are to improve the convergence properties with respect to domain of 
convergence and rate of convergence. Some results are obtained in [T5] employing experimen- 
tal approach; we supplement them with a theoretical study of convergence of the transformed 
series. 

We reconsider the derivation of (11. 3p . trying the ansatz 

W = \nz - \n(p + \nz) + u . (4.1) 
Substituting into the defining equation We w = z, we obtain 

ze u 



In z — \n(p + In z) + u 



p + In z 



From this, it is clear that if we define 



1 p + ln(p + lnz) 

a = — and r = - , (4.2) 

p + mz p + mz 

then we recover the equation (11.41) originally given by de Bruijn for u and leading to the 
series fll.3p . Thus the fundamental relation f ll.4p is invariant with respect to p, with only the 
definitions of a and r being changed. This remarkable property is due to a similarity property 
of the solution (11. ip of the original transcendental equation y a e y = x and explained by the 
following theorem. 

Theorem 22. An untrans formed series solution of equation y a e y = x and the corresponding 
transformed series for the Lambert W function, defined by (I4.2p . are connected by relations 

x = (az) a , a = e p . (4.3) 

Proof. The solution (jl.ip of equation y a e y = x possesses a similarity property with respect to 
parameter a > in the sense [17] 

/ T l/a\ / r l/a\ 

® a {x) = a<S>A X —\=aW^—\ . (4.4) 
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It follows from (EQ and (@3D that 



W 



,1/a 



a 



1 - T 



a 



U 



(4.5) 



where a and r are defined by (11.21) . The right-hand side of (14.51) does not include a explicitly. On 
the other hand, a is included in the left-hand side through a combination z = x l / a /a. Therefore, 
the fundamental relation ( 11. 4p will retain if we change variable x = (az) a . Substituting this 
formula into ( 11. 2p and introducing parameter p = lna we obtain exactly equations (14. 2p and 
the theorem follows. □ 



Thus introducing the invariant parameter p generates an infinite one-parameter family of 
series formed by replacement of variables r and a in the original series by expressions ( 14. 2p . 
Similar series for W are associated with an invariance observed in [17] and studied in [TO] . 

We now consider the properties of the transformations (14. 2 p for zeK. We shall start with 
pGK and later consider briefly one complex value of p. Both o and r are singular at z s = e~ p , 
with the special case p = recovering the previous observations regarding the singularities 
at z — 1. We note a is monotonically decreasing on z > z s . For r, we have t{zq) = at 
zq = exp(z s —p), with r positive for larger z and negative for smaller. Also we note that r has 
a maximum at z = exp(ez s — p). In Figure [3j we plot a and r, defined by (14.21) . for different 
values of p. We see that for all z > z s , a decreases with increasing p, but r increases. In view 
of the form of the double sums above it is not obvious whether convergence is increased or 
decreased as a result of these opposed changes. This is what we wish to investigate here. 



- p-0 
-P-2 



T(/,|l) 



2 



4 6 



Figure 3: Dependence a and r on z for different values of parameter p. 



4.2 Domain of convergence 

We wish to investigate first the domains of z G M for which the series (11. 3B and (11.51) converge, 
and how the domains vary with p. We begin with theoretical results. 

For p = the domains of convergence are known from theorems [2] and [TTJ Specifically, the 
series (I1.3P converges for z > e and the series (II. 5p converges for z > zq = 1.004458... (see 
Corollary [141) . For arbitrary real p the following statement can be proved for the series (11.31) . 
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Theorem 23. The domain of convergence of the transformed series (jl.3p is defined by equations 



9fcW_i (-e^ 1 (p + \nz)) > p - 1 and z > e~ p , (4.6) 

which is equivalent to 

z>z p = { ' y ~ 4.7 

I e~ P+r?oCSCJ?0 p > 

where r] G (0, 7r) zs i/ie root of equation rj cot r] = 1 — p. 

Proof. The proof of the theorem is similar to that of Theorem [2] and based on an application of 
Theorem [T] to the transformed series (II. 3p . In particular, substituting the expressions ( 14. 2 p in 
( 12. ip we obtain in the real case, i.e. under assumption p+ln z > 0, the inequality ( 14. 6 p . Applying 
Lemma [1] to the latter we get (14. 7p . where z p > e~ p , which justifies the above assumption and 
the theorem follows. □ 

Remark 24. The formula (14.71) can be obtained by substituting equations (14. 3p in (I2.14p . 

However, in contrast to x a defined by (I2.14p . z p is monotone decreasing with p. Thus, the 
larger p, the wider domain of convergence of the transformed series (jl.3p (cf. Remark^. 



Remark 25. In the formula (14. 7p . when p > but p ^ 1, we can also write z v = e 



-p+(l-p) sec?7o 



Remark 26. The convergence condition (14. 6p can be extended to the case of complex z similar 
to the extension of the condition (I2.16P for the untransformed series (II. 3p by Theorem [3 

To find out the domain of convergence of the transformed series (jl.5p one should substitute 
(14. 2 p in (13. ip and solve the obtained equation for z as a function of p. As a result, we obtain 
formulae to compute z p similar to those stated in Theorem [13j We also can find a very good 
approximation for z p for p < 1 using equation x a « x±, stated in Remark [T6l and Theorem 1221 
Specifically, substituting (14. 3p in this equation we obtain 

z p « e~ p (z o y~ P . (4.8) 

The accurate and approximate values of z p are depicted in Figure H] by solid line and circles 
respectively together with curve (14 .7p (dashed line). It follows from FigureHJas well as from (14.71) 
and (14. 8p that with increase of parameter p the domain of convergence of the transformed series 
monotonely extends. To illustrate and qualitatively verify this result we design an appropriate 
numerical procedure. The method is simply to compute the partial sum of a series to a high 
number of terms, using extended floating-point precision as necessary, and then to plot the 
ratio of the partial sum to the exact value (the exact value is obtained using a built-in Maple 
function LambertW(k,x), where a method different from series summation is used). The edge 
of the domain of convergence is then signaled by rapid oscillations and by marked deviations 
from the desired ratio of 1. (To make a graph be readable we depict only the relevant part of 
each curve.) 

For the series (ll.3p we have plotted in Figure the partial sum to 40 terms for different 
values of p. For p = 0, we see a nice illustration of Theorem [21 with the partial sum becoming 
unstable in the vicinity of z — e. For positive p, we see the domain of convergence increased 
and for negative p it is decreased, in accordance with Theorem [23j Similar effects can be 
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Figure 4: Behavior of boundary of convergence domain as a function of p for series (11. 3p (dashed 
line) and (11.51) (solid line) in real case. 

seen for (11.51) . we plot in Figure [6] the partial sums for 40 terms as p varies. The domain of 
convergence for each p is clearly seen, and confirms that the point of divergence moves to larger 
z for decreasing p and to the left for increasing p. For p = this point is very close to 1, which 
sharp demonstrates the result in Theorem [TT1 

We can summarize the above findings by noting that series (11.51) has a wider domain of 
convergence, and a better behaviour with p, while the domain of convergence for series (11.31) 
becomes worse in that order. 

The fact that the domain of convergence of the transformed series is extending while the 
parameter p is increasing can also be found in the complex case based on the results of theorems 
[20] and [231 To make certain of this it is sufficient for the series (II. 5p . to substitute expressions 
(14. 2p (with z G C) into equation (13.11) and for the series (11.31) . to consult Remark [261 The 
results are presented for p = — 1,-1/2, 0, 1/2 and 1 in Figure [7] and Figure [8] for the series (11.31) 
and (ll.5p respectively where the curves for p = are the same as in Figure [2] and the points of 
intersection of the curves with the positive real axis correspond to points on the curves depicted 
in Figure HI 

4.3 Rate of convergence 

By rate of convergence, we are referring to the accuracy obtained by partial sums of a series. 
Given two series, each summed to N terms, the series giving on average a closer approximation 
to the converged value is said to converge more quickly. The qualification 'on average' is needed 
because it will be seen in the plots below that the error regarded as a function of z can show 
fine structure which confuses the search for a general trend. Further, the comparison of rate 
of convergence between different series can vary with z and p. For some ranges of z, one series 
will be best, while for other ranges of z a different series will be best. Although one series 
may converge on a wider domain than another, there is no guarantee that the same series 
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Figure 5: For series (II .3p . the ratio W^°\z,p)/W(z) as functions of z for p = —1/2, 0, 1. 

will converge more quickly on the part of the domain they have in common. The practical 
application of these series is to obtain rapid estimates for W using a small number of terms, 
and for this the quickest convergence is best, but this will be dependent on the domain of z. 

The previous section showed that positive values of the parameter p extend the domain of 
convergence of the series, but its effect on rate of convergence is different. Figures 191 and ITUl show 
the dependence on z of the accuracy of computations of the series (I1.3P and (jl.5p respectively 
with iV = 10 for p = — 1, —1/2,0 and 1. One can see that the behaviour of the accuracy is 
non-monotone with respect to both z and p although some particular conclusions can be made. 
For example, one can observe that for the series (ll.3p at least for z < 30 within the common 
domain of convergence the accuracy for p = —1/2, and 1 is higher than for p = — 1. For the 
series (jl.5p an increase of positive values of p reduces a rate of convergence within the common 
domain of convergence i.e. for z > 1.5. However, at the same time for z > 11 computations 
with p — — 1 are more accurate than those with positive p and for 5 < z < 18 the highest 
accuracy occurs when p = —1/2. 

The next two figures [TT] and d2] display the dependence of convergence properties of the 
series (11 .3p and (II. 5p respectively on parameter p for different numbers of terms N = 10, 20 and 
40. Again, the curves in these figures confirm that the accuracy strongly depends on parameter 
p and is non-monotone and show that on the whole an increase of the number of terms improves 
the accuracy. It is also interesting that there exists a value of p for which the accuracy at the 
given point is maximum; this value depends very slightly on N and approximately is p ~ —0.75 
in Figure [TT1 and p ~ —0.5 in Figure [121 

The explained behaviour of the accuracy depending on parameter p shows that introducing 
parameter p in the series can result in significant changes in accuracy. The pointed out non- 
monotone effects of parameter p on a rate of convergence can be due to the aforementioned 
non-monotone behaviour of r. 
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Figure 6: For series ( II. 5p . the ratio W^°'(z,p)/W(z) as functions of z for p = —1,0, 1. Com- 
pared with Figure [5j this shows convergence down to smaller z. 
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Figure 7: Domains of convergence of series (11.31) in complex z-plane for p = 
-1,-1/2,0,1/2 and 1. 

5 Branch —1 and complex p 

The above discussion has considered only real values for the parameter p. We briefly shift our 
consideration to complex p and to branch —1. For z in the domain — 1/e < z < function 
W-i(z) takes real values in the range [—1, — oo). The general asymptotic expansion takes the 
form [8] 

W- X (z) = ln(» - 2ni - ln(ln(» - 2ni) + u , (5.1) 

where u denotes the remaining terms similar to (II. ip . This will clearly be very inefficient for 
z G [—1/e, 0) because each term in the series will be complex, and yet the series must sum to a 
real number. If, however, we utilize the parameter p, we can improve convergence enormously. 
We again adopt the ansatz used above to write 



W-i(z) = [ln_i z + p] — [p + ln(p + ln_i z)] 



p + ln(p + ln_i z) 
p + ln_! z 



+ v 



(5.2) 
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Figure 8: Domains of convergence of series (11.51) in complex z-plane for p = 
-1,-1/2,0,1/2 and 1. 




-9 J 



Figure 9: For series (11. 3ft with N — 10, changes in accuracy in z for p — — 1,-1/2, and 1. 

where v stands for the remaining series which will not be pursued here. By setting p = in, 
we can rewrite [ln_i z + m\ as ln(— z). When p = 0, (15. 2p is equivalent to (15. ip . A numerical 
comparison of partial sums can be used to show the improvement. Specifically, we compare the 
following first terms in ( 15. II) and ( 15. 21) 

m , , s . . . , ln(ln(,2) —2ni) 
Wr{ = ln z - 2vri - In In (z) - 2m) + ) \ J , 5.3 

m(z) — ZTTZ 

W_! = ln(-z)-ln(-ln(-z)) + ln( ~) n( ~^ ) . (5.4) 

ln(— z) 

The results are shown in Table [1] and graphically in Figure [131 One can see that both trans- 
formed (p = 27r) and untransformed (p = 0) series have an error that increases as z — > — 1/e. 
Although the maximum error occurs at z = —1/e + 0, the transformed series is exactly correct 
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H W(z) 




Figure 10: For series ( 11. 5p with N = 10, changes in accuracy in z for p = — 1,-1/2, and 1. 




1 1 ■ 



/VY v ) 



N = 


10 


N = 


20 


N = 


40 



Hi 



^ (Ar) (18,p) 



Figure 11: For series (jl.3p . the accuracy as a function of p at fixed point 2; = 18 for iV = 10, 20 
and 40. 
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Figure 12: For series (jl.5p . the accuracy as a function of p at fixed point z = 9 for N = 10, 20 
and 40. 



z 


W_i(z) 






(*) 


-0.01 


-6.4728 


-6.4640 


-6.3210- 


0.04815i 


-0.1 


-3.5772 


-3.4988 


-3.4124 - 


0.3223i 


-0.2 


-2.5426 


-2.3810 


-2.5182 - 


0.5153? 


-0.3 


-1.7813 


-1.5438 


-2.0087- 


0.6621i 


-1/e 


-1 


-1 


-1.7597- 


0.7450i 



Table 1: Numerical comparison of series transformation with p = iix. 



at z = —1/e. This is due to a difference in the local behaviour of W-\ and approximation W-\ 
when z — > — l/e + 0, particularly, the former has a square-root singularity, while W-\ is regular 
there. We also note that the transformed series is asymptotically correct as z — > 0. 
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6 Series (11.61 ) 

6.1 Different representations 

The series development (11. 6p was obtained in [TOi [9] and represents an expansion of W(x) in 
powers of cr -1 = lax 

oo 

W{x) = U + fl «( ln X T (6- 1 ) 

71=1 

oo 

W(e t )=u + ^2a n t n , (6.2) 



or 



where t = Inx and [9] 



71=1 



n— 1 // -, 

'n — 1 



The formula (16. 3p expresses the expansion coefficients a n in terms of the second-order Eu- 
lerian numbers [Hi [10] . We now show that these coefficients can also be represented through 
the unsigned associated Stirling numbers of the first kind d(m, k) given by [7] 

[ln(l +v)- v] k = £;! V (-l) m + k d(m, fc)H_ (6.4) 
^— ' m! 

m=2fc 

and the 2-associated Stirling numbers of the second kind used in the series (II. 5J) . 
Both representations can be obtained on the basis of a relation 



W(e*) + lnW(e*) = t (6.5) 

and the Lagrange Inversion Theorem [5J. To apply this theorem it is convenient to introduce a 
function that is zero at t — 0. We consider function 

v = w (t) = W(e*)/u - 1 (6.6) 

and write (16. 5p as 

t = (Jo v + ln(l + v). 
Then by the Lagrange Inversion Theorem we obtain 



oo 

V = 

n=l 



tl^^^r ( ,r> 

n — 1 >■ * 



where the operator [v p ] represents the coefficient of v p in a series expansion in v. Comparing 
(l6.6p .f l6T2|) and (16. 7p leads to a formula for the coefficients a n , which after applying the binomial 
theorem becomes 



^0 r n-ll Vr-1^ fc ~ 1 + k \ [ ln ( X + ^) ~ v \ 

l + Wb)" 1 J ^ J V n- 1 J ^(l + ^ ) fc 



*7l 

n,_ 

fc=0 
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or by (J£3D 

n! ^ (l+w„)»+* 1 ' 

If instead of function (16. 6p to take 

h = h{t) = W (e 4 ) - quo - t (6.9) 
and apply the Lagrange Inversion Theorem to invert a relation 

t = u (e~ h -l)-h 
coming from (16. 5p . then we find in a similar way 

T lzl j_ h _ n (_i)*+i £i; * 



1 ^ fn + fc - 11 



(6.10) 

1 I , , \n+k K ' 

Finally, one more representation for the coefficients a n can be found in the following way. 
Let us consider a function 

^ = ^{t) = W{e t )-t (6.11) 

which is a simplified version of functions (16 .6p and (I6.9P : now one does not need to provide the 
zero function value at t — and here ip(0) = ojq. Then it follows from (I6.5P that 

t = e~ i '-ip. (6.12) 

This equation can also be obtained from the fundamental relation (11.41) by transformation 
u = ip + lnt, a = l/t,r = hit/t which follow from (II. 2p . 

Differentiating (16.121) in t and excluding the term from the result again using (16.121) 
result in an initial value problem for ordinary differential equation 

dip _ 1 
dt 1 + t + ip ' 
Searching a solution in the form of series 

oo 

#) = Wo + ^c/ (6.13) 

n=l 

by substituting it into the differential equation and equating coefficients of the same power in 
t one can finally find 

1 if n-l \ 

Cl = -— , C n = — r- I (n-lJCn-i + YjfcCfcCn-fc 1,71 = 2,3,4, ... (6.14a) 

1 + wo ^(1 + w ) y ^ J 

At length combining (I6.13p . (l6.1ip and (16. 2p gives 

a\ = 1 + ci, a n = c n for n > 2. (6.14b) 

In practice computing the expansion coefficients in (16.11) using recurrence (I6.14ap is faster and 
takes less digits to obtain a required level of accuracy than using either of (16. 3p , (16. 8p or (I6.10p 
which, however, being different representations of the same expansion coefficients, lead to some 
combinatorial relations considered in Section [71 
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6.2 Convergence properties 



The expansion (11.61) in fact represents a series of the Wright u function [TTJl [9] cj(z) = W)c( z ) (e z ), 
where )C(z) = \{^z — ir)/(2ir)~\ is the unwinding number of z. The Wright u function was 
introduced by Corless and Jeffrey [2] and studied in [271 M- When z ^ £ ± m for £ < — 1, 
co = uj(z) satisfies equation f(z,co) = where f(z,u) = u + lnw — z (cf. (16. 5p ). Applying 
the same approach as in Section 13.11 to this equation one can obtain the same results as in 
[27] [TO], [9] . Specifically, the nearest to the origin singularities are [10] 

Z\ — — 1 — m and Z2 = —1 + iir. (6.15) 

Note that they are connected with the singularities (13 .6p of function u = u a (r), defined by (11.41) 
or (l3T27j) . through function fl3T28|) (cf. Remark [HD 

z l = s(a,ri iy ) and z 2 = s(cr,ri 0) ). 



Thus the radius of convergence is a/1 + 7r 2 [10] and the domain of convergence is defined by 

\a\ > u (6.16) 
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The estimation of uj in the vicinity of the singularities (I6.15P is j27J |9] 



z 



lu ~ — 1 — sga(^szk)V2zk yl j as z — > z^, (fc = 1, 2) 

As in Section 13. 3[ using the Darboux's theorem one can find the asymptotic expression for the 
expansion coefficients in (16.1 ft 



. 2n- 1 

. — lrsin arctan7r 

2 v ; V 2 

712(1 + 7T 2 ) - ~ 



as n — > oo. (6-17) 



Thus, as in case of the series (II .51) for positive a (see (13.201) ). the expansion coefficients in the 
series (11.61) disclose decaying oscillations in their behavior for large n. 

In real case inequality (I6.16P read as exp(— Vl + tt 2 ) < x < exp(\/l + % 2 ). Thus from 
the point of view of the domain of convergence the series (I1.6P takes an intermediate place 
between the expansion of W(x) at the origin [8] W(x) = X]^Li( — n) n ~ 1 x n /nl^ which is valid for 
— e _1 < x < e _1 , and the series (II. 5ft which is valid for x > x\ = 1.004458... (see Corollary [141) . 
These three expansions put together cover the entire region of definition of W(x). 



7 Combinatorial consequences 

In this section we collect some combinatorial consequences resulting from the above obtained 
expressions for the expansion coefficients. 

Equating the right-hand sides of (I3.18P and (I3.33j) one can find 
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± Q (1 + X) n ~ k ~ 1 X k = ± {" + *}^'-, (7.1a) 

where summation in the right-hand side starts with one as {q} >2 = [11]. Setting /x = A/(l — A) 
in (17. lap we also find 

The identities (17. ip were obtained by L.M. Smiley in a different way in [23], where notation 
{{}} was used instead of {} >2 , and referred to as the Carlitz-Riordan identities [25]. Applying 
the binomial theorem to (17. ip leads to a pair of identities expressing the 2-associated Stirling 
numbers of the second kind through the second-order Eulerian numbers and inversely 

n + q\ fn — p — l\ // n 

p=0 



q J> 2 \g-p- 1 J \p 



n\ ^ . fn—p — l\ (n + p + 1 



p=0 



„ y"(-i)^ p i 

q V q-p ) I P + 1 



>2 



Some estimates can also be obtained by comparing the found asymptotic formulas (I3.20p 
and (13.2 ip with the explicit expressions for the expansion coefficients in (ll.5p . For example, 
taking estimate (I3.2ip and the expansion coefficients in (II. 5p at a = — 2 we obtain 

Efp + m — 11 (m — 1)! ,_ . 

< y > r as m^oo 7.2 

p=1 I P i>2 2 v /^(21n2-l) m ^ 

where the term with p = is skipped (cf. (I7.1al) . This result is consistent with the formula 
given in Ex. 10 (7), p. 224]. 

Another consequence is obtained by taking the expansion coefficients in (I1.5P at o = 
together with (I3T22D 

m_1 C i 1 

_ I T) -I- m — I ] 

(m-1)! (7.3) 



v- ( _ 1)P + m -ifP + m - 11 



Further, comparing (16.31) . (16. 8p and (I6.10p between one another we come to the following 
three identities 

(l + ^-^W k // y } ^ (l + u;o) fc 1 ' 

^//n-l\ .vfc+i , fc+ i_vi Jn + fc-ll (-l)*w§ 
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n-1 




d(n + k — 1, k) 



(n + k-U (-l) k ool 
\ k ) >2 (l + u ) k 




(7.4c) 



k=0 



Finally, combining either of ( I6.3p . ( 16. 8 p or ( 16.1 Op with (I6.17P gives an asymptotic expression for 
the sum involved there. 

Thus studying expansion series for W function we, on the way, derived the Carlitz-Riordan 
identities (17. ip as well as found a formula for an alternating sum of 2-associated Stirling numbers 
of the second kind (17. 3p and confirmed the asymptotic formula ( 17. 2ft for summation of the same 
numbers without the alternating factor. We also found formulas ( 17. 4p where the Omega constant 
Uq plays a role of a 'magic' number which connects sums involving the second-order Eulerian 
numbers, the associated Stirling numbers of the first kind and the 2-associated Stirling numbers 
of the second kind. 

8 Concluding remarks 

We ascertained the domain of convergence of the series (11.31) and (II ,5p in real and complex cases 
and found that the series ( II .5p has a much wider domain of convergence than that of the series 
( II. 3p in both cases and provided an analysis of this fact in real case. We found asymptotic 
expressions for the expansion coefficients and obtained a representation of the series (ll.5p in 
terms of the second-order Eulerian numbers. 

We considered an invariant transformation defined by the parameter p and applied it to 
the above series. We studied an effect of parameter p on the convergence properties of the 
transformed series theoretically and numerically and found that an increase of p results in 
an extension of the domain of convergence of the series. Thus the series obtained under the 
transformation with positive values of p have a wider domain of convergence than the original 
series does. However, at the same time a rate of convergence can be found to be reduced when 
the parameter p increases. Therefore in such a case within the common domain of convergence 
of the series with different positive values of p the series with the minimum value of p would 
be the most effective. 

We also considered the well-known expansion of W(x) in powers of lnx and gave an asymp- 
totic estimate for the expansion coefficients. We found three more forms for a representation of 
the expansion coefficients of the series in terms of the associated Stirling numbers of the first 
kind ( 16. 8p . the 2-associated Stirling subset numbers (I6.10p and iterative formulas (I6.14p . Finally 
we presented some combinatorial consequences, including the Carlitz-Riordan identities, which 
result from the found different forms of the expansion coefficients of the above series. 
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